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We investigate the distribution of instanton sizes in the framework of a simplified model for ensembles of instan- 
tons. This model takes into account the non-diluteness of instantons. The infrared problem for the integration 
over instanton sizes is dealt with in a self-consistent manner by approximating instanton interactions by a repul- 
sive hard core potential. This leads to a dynamical suppression of large instantons. The characteristic features 
of the instanton size distribution are studied by means of analytic and Monte Carlo methods. We find a power 
law behaviour for small sizes, consistent with the semi-classical results. At large instanton sizes the distribution 
decays exponentially. The results are compared with those from lattice simulations. 
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1. INTRODUCTION 

Instantons are field configurations of non- 
abelian SU(iV) gauge theories, which lead to non- 
pcrturbative effects. In recent years they have 
been studied in lattice gauge theories by means 
of Monte Carlo calculations by different groups 

Hi- 

In the dilute gas approximation the logarithm 
of the partition function contains an integral over 
instanton sizes, Jdp p~ 5 (pA) b , where b = 11AT/3 
and A is the scale parameter. The integrand in- 
creases with p, leading to an infrared divergence. 
This is an artifact of using the semiclassical ap- 
proximation after it has become invalid, i.e. for 
p > 1/A. If the semiclassical approximation is 
meaningful at all, a solution of this problem in 
the context of the full instanton ensemble is re- 
quired. 

The simplest way is to cut the integrations off 
at some ad-hoc value p c . But the dominant con- 
tribution comes from large pj near the cut-off 
where the assumption of diluteness fails. More- 
over, the introduction of an ad-hoc cut-off leads 
to inconsistencies with the renormalization group 
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In order to solve the problem it has been pro- 
posed that instanton sizes are cut off in a dynam- 
ical way [||]6|. The cut-off should originate from 
configurations where instantons start to overlap. 
The interaction is expected to suppress overlap- 
ping instantons and to result in a self-consistent 
cut-off. 

In connection with the dynamical cut-off the 
distribution of instanton sizes is of central impor- 
tance. For small sizes the distribution is predicted 
to be 



n(p) 



„b-5 



(1) 



by the dilute gas approximation. For large sizes 
p, where the dynamical cut-off is in effect, not 
much is known about the distribution. There are 
arguments in favour of a suppression like 



n{p) ~ exp(— cpP) with p = 2 . 



(2) 



*Talk presented by G. Munster 



We have investigated Pi the distribution of in- 
stanton sizes in a model [JfJ where the instanton 
interactions are approximated by a repulsive hard 
core potential. The radius of an instanton core 
varies proportional to the size pj of the instan- 
ton. Although this approximation appears to be 
crude, the general features of the instanton en- 
semble with a dynamical cut-off are present. 
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2. SIZE DISTRIBUTION 

In d = 1 dimensions the distribution of instan- 
ton sizes is exactly given by 

C 



n(p) = j p 



b-2 e -cp_ 



(3) 



where we recognize an exponential suppression of 
large instanton sizes. 

In higher dimensions, d > 1, we obtained ap- 
proximate expressions by means of a van der 
Waals type approximation 

Cd 

n{ P ) = T p"- 4 - 1 exp(-c/) . (4) 

For small p it grows powerlikc with the semiclas- 
sical exponent a = b — d — 1. The value of the 
exponent p = d in the exponential decay at large 
p should be considered with reservations, because 
the saddle point approximations are of uncertain 
quality there. 

We have studied the size distribution also by 
means of grand canonical Monte Carlo simula- 
tions of the simplified instanton gas model. In 
the case of d = 1 dimensions the Monte Carlo 
data agree very well with the available exact re- 
suit (g). 

In the more interesting case of four space-time 
dimensions (d = 4) we consider a = 7/3, which is 
the value for SU(2) gauge theory. The resulting 
size distribution shows the expected behaviour. 
For small instanton radii a power law with expo- 
nent a can be confirmed. In order to study the 
behaviour of n(p) for large p we considered the 
ratio F(p) = n(p)/p a and tried fits of the form 
Ffit(p) — aexp(— cp p ). In agreement with the 
theoretical results they showed that c depends on 
a, while p is nearly independent of it. In Fig. |l| 
the result of a fit in the interval [0, 2.25] is shown. 

The main interest is in the exponent p. Wc 
find a « 0.89, and the fit leads to c = 3.3 ± 0.2 
and p = 1.9 ± 0.2. For a = 6, the SU(3) case, 
the results for p are the same within the present 
errors. 

In recent years much effort has been devoted 
to lattice Monte Carlo calculations of properties 
of the instanton ensemble, and some quantitative 
statements have been given. For small p, lattice 
calculations appear to support the power law (g). 
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Figure 1. F(p) in d = 4 dimensions from a Monte 
Carlo simulation in comparison with the fit. 

For the large-p distribution, de Forcrand et al. 
predict an exponential decrease with p = 3 ± 1 
from their SU(2) lattice data 0. In contrast to 
this, Smith and Teper conclude form their SU (3) 
simulations a decay according to p~^ with £ w 
10...12|. 

To conclude, fits to our numerical Monte Carlo 
results suggest a behaviour like 



n(p) 



p a exp(-cp 2 ) . 



(5) 



The results indicate that our simplified model 
reproduces the main features of instanton ensem- 
bles with a dynamical infrared cut-off. 
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